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ANALYTICAL METHOD FOR STEADY STATE HEAT TRANSFER IN 
TWO-DIMENSIONAL POROUS MEDIA 
by Robert Siegel and Marvin E. Goldstein 
Lewis Research Center 

SUMMARY 

A general technique has been devised for obtaining exact solutions for the heat 
transfer behavior of a two-dimensional porous cooled medium. Fluid flows through the 
porous medium from a reservoir at constant pressure and temperature to a second res- 
ervoir at a lower pressure. For the type of flow involved, the surfaces of the porous 
region that are each at constant pressure are boundaries of constant velocity potential. 
This fact is used to map the porous region into a strip bounded by parallel potential lines 
in a complex potential plane. The energy equation, derived by assuming the local ma- 
trix and fluid temperatures are equal, is transformed into a separable equation when its 
independent variables are changed to the coordinates of the potential plane. This allows 
the general solution for the temperature distribution to be found in the potential plane. 
The solution is then mapped into the physical plane to yield the heat transfer character- 
istics of the porous region. An example problem of a porous wall having a step in 
thickness and a specified surface temperature or heat flux is worked out in detail. 


INTRODUCTION 

There are a number of flow and heat transfer applications that involve porous media. 
These include drainage of water in soil, the use of packed bed heaters, drying of mate- 
rials, and transpiration cooling. In the last of these, with which this report is con- 
cerned, a coolant is made to flow through the porous material in order to protect this 
material when it is exposed to a high temperature environment. Possible applications 
for transpiration cooled materials are in rocket nozzles, gas turbine blades, and por- 
tions of vehicle surfaces in high speed flight. 

Some of the early analyses of the heat transfer characteristics of porous cooled 
materials carried out around 1950 have been given in the heat conduction text by 


Schneider (ref. 1). Some of the more recent papers are reference 2, which contains a 
long list of references, and reference 3. As can be seen from these references, the 
analyses of the heat flow within a porous cooled material have all been for one- 
dimensional geometries such as for a plane slab or radial flow through a tube wall. The 
objective of the present report is to present a method for obtaining analytical solutions 
in two-dimensional geometries. 

It is evident that the configurations used for such devices as porous cooled turbine 
blades or rocket nozzles will deviate from a one- dimensional geometry. In some in- 
stances it may be possible to make a locally one-dimensional assumption in analyzing 
these geometries. However, the limits for making such an assumption can only be eval- 
uated by examining solutions in more than one dimension. When the cross section of the 
material deviates considerably from a one-dimensional shape, the fluid flow in the por- 
ous matrix can follow complex paths, and two- or three-dimensional solutions are re- 
quired. General two-dimensional solutions will be obtained here for either an arbitrary 
temperature variation or an arbitrary heat flux variation on the surface of the porous 
cooled medium. The relations that are found between surface temperature and heat flux 
would enable the solution for the heat transfer in the porous material to be coupled to the 
solutions to an external heat transfer problem such as a boundary layer flow. 

The present analysis will utilize the fact that in many instances within the pores of 
a- porous material the fluid and matrix material are in good thermal communication. As 
a result the local fluid temperature is the same as the local temperature of the matrix 
material. A single energy equation can then be written which is composed of two terms, 
one of which represents the energy carried by the flowing coolant and the other being the 
energy flow due to heat conduction in the matrix material. 

The velocity that appears in the energy equation is a function of position in the 
medium, and for the slow viscous flow encountered in the pores of many porous media, 
the velocity is proportional to the local pressure gradient (Darcy's law). By suitably 
defining a dimensionless velocity and pressure, the dimensionless pressure is fixed at 
zero on the coolant reservoir side of the porous medium, and at unity on the high tem- 
perature side of the medium. Since the dimensionless velocity is proportional to the 
gradient of the dimensionless pressure, and since in addition, the velocity satisfies the 
continuity equation, the dimensionless pressure can be regarded as a velocity potential. 
Thus the region occupied by the porous material can be mapped into the region of a com- 
plex potential plane lying between two parallel potential lines at zero and unity which 
correspond to each of the constant pressure surfaces. 

The method given here is based on transforming the energy equation so that its inde- 
pendent variables become the coordinates of the potential plane. In this plane the trans- 
formed energy equation is separable and the geometry is in the simple form of a unit 
strip. These facts make it possible to obtain a general solution to the energy equation. 
After the general solution has been obtained in the potential plane, conformal mapping is 
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used to relate it to the coordinates of the physical plane. 

The method described can be carried out for a variety of geometries and boundary 
conditions. The general solution is obtained here as a typical case for a two- 
dimensional porous region that is long in one of the coordinate directions and has a vari- 
able thickness in the other direction. One surface of the medium can have either an 
arbitrary specified temperature or heat flux. 

To illustrate the application of the general solution, the specific situation of a step 
porous wall (i. e. , a wall that changes in a step function fashion from one constant thick- 
ness to another) is considered. The solution is carried out to yield surface heat fluxes 
or temperatures corresponding to cases where one surface of the medium is either 
maintained at a specified uniform temperature, has a step function temperature varia- 
tion, or has a specified uniform heat flux. 


SYMBOLS 


^ 2 * ^ 4 ’ 

/V 

E 


E 


F 

G 


H 


h r 

i, 3 

k m 

L s 

Z s 

M 


N 

n 


P 

q 


ratio of thick to thin dimensions of step porous wall 

integration constants 

specific heat of fluid 

transform defined by eq. (44) 

function defined in eq. (65b) 

function in specified temperature distribution eq. (6a) 

function in specified heat flux distribution eq. (6b) 

transform defined by eq. (39) 

reference length in porous material 

unit vectors in X- and Y-directions 

thermal conductivity of porous region 

dimensionless coordinate along boundary S, Ij h 

o r 

coordinate along boundary s such that l - 0 at x = 0 

s 

heat flux parameter, (qg - qj)/q^ 
temperature parameter, (tg - tj)/(t^ - t^) 
outward normal vector 
pressure 
heat flux 
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T 

t 
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u 

V 

v 


a 
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X,Y 


x,y 

z 

fi 

A 

6 

s 

v 

e 

K 

X 

M 

I 

P 

$ 

X, a) 
* 


bounding surfaces of porous region in dimensionless coordinate system 

bounding surfaces of porous region 

boundary surface of arbitrary volume -a- 

dimensionless temperature defined in eqs. (10) and (11) 

temperature 

it 

dimensionless velocity u 

K Po - Ps 

velocity 

dimensionless velocity in Y direction 

velocity in y direction 

arbitrary volume 

complex potential W = + icp 

dimensionless coordinates, — and 

h r h r 

rectangular coordinates 
dimensionless physical plane, X + iY 
separation constant in solution for 9 
quantity defined in eq. (11) 

Dirac delta function 

intermediate mapping plane, £ = £ + irj 

imaginary part of ? 

dependent variable defined by eq. (29) 

permeability of porous material 

P c D «(p - p ) 
parameter, — s~ - 

2k m 
m 

fluid viscosity 
real part of £ 
fluid density 

function of <p in solution for 9 
potential, imaginary part of W, (p Q - p)/ (p Q - P s ) 
dummy variables of integration 
function of \f/ in solution for 6 
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\f/ real part of W 

V dimensionless gradient, eq. (15) 

Superscript: 

at location of step in plate thickness 
Subscripts: 

o at the boundary s Q or S Q 

s at the boundary S 

°° coolant reservoir 

1,2 values at large negative and large positive X (or x) on surface S (or s) 


DERIVATION OF BASIC EQUATIONS 


Let (A be any volume within a porous medium with effective thermal conductivity 
k m (based on the entire cross sectional area of the porous material) and permeability k 
such that u is large with respect to the pore size and let i. denote the surface of o-. 
Suppose that there is a fluid with constant density p, constant heat capacity Cp and 
constant viscosity p which is flowing through ia. Assume that the thermal conductivity 
of the fluid is very small compared with k m and that the pore size is so small that the 
fluid obeys Darcy’s law. Let u denote the local Darcy velocity of the fluid (this is the 
velocity obtained by dividing the volume flow by the entire cross section of the porous 
material rather than the open area). If the thermal communication between the fluid and 
the matrix is sufficiently good, the local fluid temperature will be approximately equal 
to the matrix temperature in the immediate vicinity. We denote this common tempera- 
ture by t. Finally, suppose that a steady state situation exists within ia. Then an 
overall energy balance applied to o- shows that 


J (pCptu - k m Vt) • d? = 0 


or applying the Divergence theorem 


J [pC p V • (tu) - k m V 2 t]dP = 0 


If the changes in t and u are very small in distances on the order of the pore size, 
we can conclude from this (since o- is arbitrary) that 


5 



II 


k m v2t " pC p V ' (S) = ° (1) 

Now the equation of continuity for the fluid is 

V • u = 0 (2) 

This shows that equation (1) can be written as 

k m v2t ' pC p" ' Vt = 0 (3) 

Darcy's law for the fluid velocity is 

u=-iVp (4) 

where p is the local pressure of the fluid within the medium. Equations (2) to (4) are 
the three-dimensional generalizations of the classical one- dimensional porous cooling 
equations (see ref. 1). 

GENERAL ANALYSIS OF TWO-DIMENSIONAL POROUS COOLED WALLS 

The steps in the analysis involve (a) establishing the boundary conditions associated 
with the porous wall, (b) formulating these boundary conditions and the associated energy 
equation in dimensionless form, and (c) transforming the entire boundary- value problem 
into the potential plane and finding a general solution. The general solution applies to 
any arbitrary temperature or heat flux distribution along one surface of the porous mate- 
rial. Conformal mapping is used to relate the potential plane solution to the physical 
plane. 


Boundary Conditions on Porous Region 

Now consider the two-dimensional porous wall shown in figure 1. The lower surface 
of the wall whose unit outward-drawn normal is n Q is denoted by s , and the upper sur- 
face of the wall whose unit outward-drawn normal is n is denoted by s. We suppose 

s 

that no changes occur in the direction perpendicular to the x, y-plane. Below the wall 
there is a reservoir which is maintained at constant pressure and temperature p Q and 
t^ , respectively. The pressure in the fluid above the wall is constant and equal to p g . 
We suppose that p Q > P s - Then the fluid flows from the reservoir through the porous 
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Figure 1. - Two-dimensional porous cooled wall. 


wall and out through the top surface. Since p Q and p g are both constant, the fluid ve- 
locity at the wall surfaces (s Q and s) will be in a direction perpendicular to these sur- 
faces. The conditions on the upper surface of the wall will be chosen so that heat flows 
by conduction through the wall from the upper surface into the reservoir below at t^ . 

We assume that the assumptions of the previous section apply within the wall and there- 
fore that the flows of heat and mass within the wall are governed by equations (2) to (4). 

Now as the fluid in the reservoir approaches the wall, its temperature will rise 
from the reservoir temperature to the wall surface temperature t Q , which is an 
unknown variable along s Q . Since the thermal conductivity of the fluid is assumed to be 
much less than the thermal conductivity of the wall, for reasonable magnitudes of the' 
fluid flow velocity the thickness of the fluid layer over which this temperature rise takes 
place is very small compared with the wall thickness. We can therefore assume that the 
fluid layer is locally one dimensional. Since the velocity is perpendicular to s Q , there 
is no flow along this thermal layer. Hence, applying an energy balance to the thermal 
layer shows that 


k mV Vt= P C p (t o -t «> )fi o- " 


P = P r 


for (x, y)es 


( 5 ) 


We shall consider two different types of boundary conditions for the upper surface 
s. First, we shall suppose that its surface temperature is specified so that it varies 
along the surface from a temperature t^ at £ g = -°° to a temperature tg at = +°° . 
The Z is the distance measured along the upper surface from the point where it inter- 

o 
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sects the y-axis. Thus in this case the boundary conditions on the upper surface of the 
wall are 


t =t Q = t, + (t 9 - t. )FI — — 


”s " "1 
P - Pc 


for (x, y)es 


(6a) 


where F is a given function which is equal to zero at[ — )= -«> and equal to 1 at 

/i s \ ^ 

I — I = +°° , and h r is a reference length. 


The second boundary condition which we shall consider is where the heat flux into 
the wall is specified along the upper surface so that it varies from a value of at 
Z g - -co to a value of q 2 at Z g = +°° . Under these conditions the boundary condition 
on the upper surface becomes 


k m S s ' vt * <ls = + (q 2 - 

' u r' r 

p = p c 


for (x, y)es 


(6b) 


where G is a given function which is equal to zero at Z g /h r = -°° and equal to 1 at 

l /h = +°° . Equations (2) to (4) together with the boundary conditions (5) and (6a) ot 
s r 

alternatively together with the boundary conditions (5) and (6b) completely determine the 
solution to the heat transfer problems within the porous wall. 


Dimensionless Form of Energy Equation and Boundary Conditions 

It is now convenient to introduce the following dimensionless quantities: 


= p c p k K - p s > 

(7) 

2k ju 


m 


N = * 1 

(8) 

r-b 

1 

rb 

8 


q 2 ~ qi 

M = ^ 1 

(9) 


*1 
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where 


X = x/h r 
Y = y/h r 


L s = V h r 



U = M. 


r u 


* (P 0 - P S ) 





>- 


J 


( 10 ) 


(tj - ) if boundary condition (6a) applies 


A = 


1l h r . 


(ID 


^. k m 


if boundary condition (6b) applies 


Upon substituting these definitions into equations (2) to (4) and boundary conditions 
(5), (6a), and (6b), we obtain 


V 2 T - 2 XU • VT = 0^ 


U = Vcp 


V • u = 0 


(12) 


n Q • VT = 2Xn o • UT 


cp = 0 


for (X, Y)eS r 


(13) 


T = 1 + NF(L f 


<P = 1 


for (X, Y)eS 


(14a) 
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or 


n s • VT = 1 + MG(L g f 

> 

<p = 1 > 


for (X,Y)eS 


where 

ax 9Y 

and the porous wall in dimensionless coordinates is shown in figure 2. 


(14b) 


(15) 



Figure Z - Dimensionless physical plane (Z-planeh 

The second equation (12) can be used to eliminate U in the other two equations (12). 
Thus 

V 2 T - 2\Vcp ■ VT = 0 (16) 

and 

V 2 <p=0 (17) 
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Since <p is constant on S Q and S, it is clear that 

s °‘Tft for(x>Y)eS ° 

and 

for (X, Y)eS 

S I *<p\ 

Using these results together with the second equation (12) in the boundary conditions (13) 
and (14b), we obtain 


and 


V<p • V T = 2AT I V cp 

cp = 0 


for (X, Y)eS r 


(18) 


Vcp • VT = | Vcp | [l + MG(L g )j 


for (X, Y)eS 


(19) 


<P = 1 


Transformation of Boundary Value Problem Into Potential Plane 


Since equation (17) shows that cp satisfies Laplace's equation, there must exist a 
harmonic function i// and an analytic function W of the complex variable 

Z = X + iY (20) 


such that 


W ^ i// + i cp (21) 

Physically the change in \js between any two points is proportional to the volume flow of 
liquid crossing any curve joining those two points. Hence \p must vary between -°° 
and +°° as X varies between - 00 and +°° . Since 

9*tiW = 0 for ZeS Q 
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and 


9mW = 1 for ZeS 


The mapping 


Z - W 

transforms the physical plane (Z-plane) shown in figure 2 into the infinite strip shown in 
figure 3 in the W-plane. The boundaries S Q and S in the physical plane are trans- 
formed into the lines <p - 0 and cp - 1, respectively, in the W-plane. The mapping 
Z — W can be found by using conformal mapping once the shape of the porous wall in the 

v 


< p = l 


0 >= 0 


Figure 3. - Complex potential plane (W-plane). 


physical plane is specified. This gives the solution to the boundary value problem for 
the flow. (An example for a specific wall geometry will be worked out subsequently. ) 

The boundary value problem for the temperature in the physical plane can be solved 
if we use the Boussinesq transform (ref. 4) to transform it into a boundary value prob- 
lem in the W-plane. Thus, the independent variables X and Y in equation (16) and its 
boundary conditions will be changed to the variables cp and , and the resulting bound- 
ary value problems will be solved in the infinite strip in the W-plane shown in figure 3. 
Once this solution, which gives T as a function of cp and ip has been found, the map- 
ping W — Z (which is completely determined once the geometry of the wall in the physi- 
cal plane is specified) can be inverted to give Z as a function of cp and \p . Thus T 
will be known parametrically (in terms of the parameters cp and i|/) as a function of X 
and Y. This will complete the solution to the problem. 
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To transform the energy equation (16) and boundary conditions into the potential 
plane, the following relations are used (ref. 5): 


V 2 T = 


/a 2 T ( 9 2 t\ 

dW 

W 2 dcp 2 / 

dZ 


( 22 ) 


and 


Vcp 


dW 

dZ 


(23) 


Notice that 


dW 

dZ 



Hence, upon taking real and imaginary parts, we find 


dcp _ 

dW 

2 ax 

ax 

dZ 

dcp 

dcp _ 

dW 

2 BY 

BY 

dZ 

dcp 


This shows that 


V cp 


/3T ax + dr dY\ 

dW 

2 

_ aT 

dW 

\aX dcp dY dcp) 

dZ 


a cp 

dZ 


(24) 


Finally, notice that since cp is constant on S, the distance L c along S is a function of 

s 

\[/ only. (This functional relation is known once the mapping W — Z, which solves the 
flow problem, is known. ) Thus 


L g (X,Y) = L s fc//) for (X, Y)eS (25) 

Now, using equations (22) to (25) in equation (16) and the boundary conditions (18), 
(14a), and (19) yields 


dxf / 2 dcp 2 d(p 


( 26 ) 
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— - 2 AT = 0 

dtp 

T = 1 + NF (L g (ip )) 


for cp — 0 
for ip - 1 


(27) 

(28a) 


or 


ST 

dip 


dZ 


dW 


[l + MG(L S (*))] 


for tp - 1 


(28b) 


Equation (26) can be simplified by introducing the new dependent variable 0 defined 


by 




(29) 


In terms of this new variable equations (26) to (28b) become 




dip 2 dip 2 


00 


xe 


dip 

9 = 1 + NF (L g (i//)) 


at (p = 0 

at tp - 1 


(30) 

(31) 
(32a) 


or 


00 _ 
0 ip 


+ xe = 


dZ 

dW 


[l + MG (L s (!//))_ 


at = 1 


(32b) 


Equation (30) together with the boundary conditions (31) and (32a) (or alternatively 
(31) and (32b)) constitute a boundary value problem for 0 in the unit strip in the W- plane 
(shown in fig. 3) which completely determines 0 (and therefore T) as a function of tp 
and x (/ . Notice that the particular relation between Z and W (which depends only upon 
the geometry of the porous wall in the physical plane) enters only through the boundary 


condition (32a) (or alternatively (32b)) since L and 


of 


dZ 


dW 


must be known as functions 


ip in order to completely determine these boundary conditions 

I dZ ' 

solve these boundary value problems for arbitrary Lg and 1 


dW 


It is possible to 
(as well as arbitrary F 


and G), and so the particular geometry and boundary conditions can be substituted into 
the general formulas once this general solution has been obtained. 
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General Solution of Boundary Value Problems in Potential Plane 


The boundary value problems posed by equation (30) and the boundary conditions 
(31) and (32a) (or alternatively the boundary conditions (31) and (32b)) can easily be 
solved by the method of separation of variables. To this end we seek a solution to equa- 
tion (30) of the form 


d{ip, (p) = *(\p)4>((p) 

Upon substituting this into equation (30), we find 

_ x 2 
<f» 


This implies that there exists a constant /3 such that 


and 


T t t 
* 



4> 

<i> 




(33) 


Hence, 


4> = 


C^e 





+ C 2 e 





(34) 


= 


C 3 e 


m 


C 4 e 


-m 


(35) 


where through are arbitrary constants of integration. 

Upon substituting equations (33) to (35) into the boundary condition (31), we find that 
this boundary condition will be automatically satisfied if 

Cj/xV - C 2 Va 2 +| 3 2 = X(C 1 + C 2 ) 


or 


Vx 2 + )3 2 - A 
A 2 + /3 2 + A 
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Using this in equation (34) gives 


$ =. 


2C, 


A 2 7? + x 


X sinh (»^7 + Vx 2 + Z3 2 cosh (Vx 2 + /3 2 X 


(v777*>)' 


Hence, the solution to both boundary value problems must be of the form 

e i/3 ^ d/3 (36) 


2 77 •'-CO 


C(/3) 


A sinh A 2 + /3 2 <pj + Va 2 + /3 2 cosh^A 2 + /3 2 (p 


where the function C of /3 can be determined so that either the boundary condition (32a) 
or the boundary condition (32b) is satisfied. 

Solution to boundary value problem for specified temperature. - First suppose that 
the boundary condition (32a) applies. Put 


H(/3) = COS) 
Then equation (36) becomes 


A sinh^V^A 2 + /3 2 ^ + Va 2 + /3 2 cosh (l/ A 2 + /3 2 ^ 


f 

»/. oo 


H(/3) 


A"_+ J3" <p)_ + r A" + /3" cosh( V A 2 + /3 2 <p 
A sinhV P A “ + (3“J + r A 2 + /3 2 cosh^A 2 + 0 2 ) 


A sinh ( Va 2 + / 3 2 <p) +_/' 2 ■ ° 2 
inh(l^ A 2 + /3 2 ) + >^' 2 ■ ° 2 




d/3 (37) 


Substituting this result into the boundary condition (32a) yields 

NF(L g (^))=-l f H(/3)e i/3 ^ d/3 

2n J - oo 


1 + 


Upon inverting this Fourier transform, we find 


HO) = f [l + NF(L g (^))]e- i ^ di// 

= f e _i ^d^ + N f F (L (^/))e" i ^^ / di// 

J- OO J- oo k 
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The first definite integral is equal to 


r e~ i(3xf/ dx(y = 2tt5(J3) 

oo 


where 6 is the Dirac delta function. Hence 

H(/3) = 2 7r 5 (/3 ) + NH(/3) 

where 


( 38 ) 


/♦oo 

H(i3)H JL F ( L s W) e ~ i/3,//d ^ 


(39) 


Substituting equation (38) into equation (37) yields 


2tt 



H(/3) 


X sinh 


t777 ^ t cosh(» / x 2 + /3 2 »>)" 

Vi 2 + 0 2 ) ♦ V ^ 2 + 0 2 co.sh(V A 2 + (3 2 ) J 


A sinh l FA + /3' 


e i/3 ^ d/3 


(40) 


Thus equation (40) with H defined in terms of the dimensionless temperature dis- 
tribution F by equation (39) is the solution to the boundary value problem. 

In this case it is also of interest to have an expression for the conduction heat flux 
q crossing into the surface S, that is, 

b 

q g = k m h g • Vt for (x,y)es 


Therefore for (X, Y)eS, by use of equations (24) and (23), 


q s h r = fi . _ V<? • VT = 

^(tj-U s |v<H 

Hence, by use of the transformation (29) and the fact that 
find 


dW dT 

dZ d<p 

cp = 1 on the boundary S we 


^s h r 

dW 

(M. + \e\ 

k m<‘l - ‘ J 

dZ 

U ) 


(41) 
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Upon substituting equations (32a) and (40) into this expression, we obtain 


% h r 


k mftl - *-0) 


dW 


dZ 




<P = 1 


2A + ANF(L (ip)) 

D 


_N 

2n 


f 


+ /3 2 


A cosh^A 2 + /3 2 ) + V^A 2 + /3 2 sinh(^A 2 + /3 2 j 


A sinh(A 2 + /3 2 ) + ^A 2 + /3 2 cos1i(Va 2 + /3 2 ) 


d/3 > 


(42) 


Solution to boundary value problem for specified heat flux . - Now suppose that the 
boundary condition (32b) applies. Put 


E(/3) = C03) 


(2 A 2 + j3 2 ) sinh(| // A 2 + £ 2 ) + 2aVa 2 + /3 2 cosh A 2 + /3 2 ) 


Then equation (36) becomes 


%/- OO 


EOS) 




XsinhlfX +/3 <£7+ ?X + /3 cosh 


2 . «2 1 . ^2 . 0 2 




(2X‘ » jS“) sinh(tx 2 + /3 2 ) + 2*VT 2 T/3 2 cosh If A‘ + [ 1 ‘ 


.2 . „2x 


■fa 


2 . n 2 


e 1/3 ^ d/3 


(43) 


Substituting this result into the boundary condition (32b) yields 


dZ 


dW 


[1 + MG(L ( 


<P= 1 


iW)l=f f 


E(j3)e i/3 ^ d/3 


Upon inverting this Fourier transform, we obtain 


EOS) = 


f 


dZ 


dW 


[1 + MG(L_ (^))]e' i ^d^/ 




(44) 


Thus equation (43) with E defined in terms of the dimensionless heat flux distribution G 

by equation (44) is the solution to this 


and the reciprocal surface velocity 


dW 


<P= 1 
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boundary value problem. In particular the upper surface temperature distribution t g 
(on S) is given by evaluating equations (29) and (43) at cp - 1 to obtain 


(t s ~ t|: p^ k m 
*^s h r 


~f 

2n / 

%/- OO 


E(/3) 


X sinh( V 

X 2 ♦ (3 2 ) 

+ 

V X 2 + /3 2 cosh 

(l'? 



(2X 2 + /3 2 )sinh ( 

(Vx 2 + 

/3 2 ) 

* + 2X^X 2 + /3 2 

COSlJ 


+ /3 2 ) 


M3* 


d/3 (45) 


SOLUTION FOR A STEP POROUS WALL 

The general solution has now been obtained for either a specified arbitrary tempera- 
ture or heat flux at one surface of the porous material. To illustrate the application of 
the solution, results will be obtained for the example of a step geometry as shown in fig- 
ure 4(a). The thickness of the thin portion of the wall has been chosen as the reference 
dimension, so that all lengths are in units of this dimension. 
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(a) Dimensionless physical plane, Z = X + iY. 
<P 

r <p = i 

c 


ip 

(b) Potential plane, W = ip + \<p. 

Figure 4. - Mapping of step porous wall into unit strip of po- 
tential plane. 
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Conformal Mapping Relations Between Z and W Planes 


In the solutions the quantity 


dW 


dZ 


<M 


appears along with the functional relation 


L g (\p) relating positions along the surface S to positions along the boundary <p = 1 in 
the potential plane. These quantities can be found from the function Z — W, which maps 
the region shown in figure 4(a) into the region shown in figure 4(b). This function can be 
obtained in terms of an intermediate variable £ by utilizing the mapping given in refer- 
ence 6. The £ plane is shown in figure 5. The origin in the physical plane corresponds 


v 

t 

|5 


1 6 1, 2 3,4 5 

£=1 <f = A 

Figure 5. - Mapping of step porous wail into intermediate £ -plane, 

C ' £ + 1 V- 


to W = 0 and £ = to. The mapping is given by 


Z =- 

[Alo^ 1 

- + A )l 

n 

Vc - 1 

/ V? - A/J 


and 


5 = 4 + irj = 



To obtain dW/dZ 


we use the relation 
<P = 1 


(46) 


' (47) 


dW _ dW dg 
dZ d? dZ 
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The term dZ/d£ is found from equation (46) to be 


dZ _ 2A A 2 - 1 

d ? 77 2 _ j )( ? 2 _ a 2 j 


Upon solving equation (47) for W, we find 


Hence, 


so that 


1 j“2 _ a 2 

W =± ln^ — 

77 ? 2 - 1 


dW _ 2? 


a 2 -i 


d ? jt ^2 _ A 2) (i *2 _ 


dW = £ = J. / e 7rW - A 2 
dZ A A \ e ?rW _ j 


, 1/2 


At (p = 1, e 77 " 7 = e 77 ^ 7 ' 77 = -e 777 ^. Therefore, 


dW 

dZ 


<P=1 


e ff ^+ 1 


(48) 


(49) 


(50) 


To obtain the relation between L and xjy along the boundary <p - 1, it is conven- 
ient to use £ as an intermediate variable. The mapping defined by equation (46) trans- 
forms the porous region into the upper right quarter £-plane as shown in figure 5. The 
boundary cp = 1 maps onto the £-axis between £ = 1 and £ = A. The relation between 
L and | is found by letting £ - £ + i rj in equation (46). Then taking the real part 
gives an expression for X as a function of £ and rj . Upon setting rj = 0 in this ex- 
pression, we obtain an equation for L as a function of The result is 

b 

L =1 (A loglil- log^ztl\ 1 < 4 < A (51) 

s jt V 4-1 A - 4/ 
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The relation between £ and j// along L g is found from equation (47) by letting 77 = 0 
and cp = 1 to obtain 




/e^+ A 2 \ 
\e'* + l / 


-00 < < 00 


(52) 


Thus equations (51) and (52) are parametric equations that relate L to xj/. 

s 


Exit Velocity From Porous Plate 


It follows from equations (12) and (15) that for this example the velocity of the fluid 
leaving the porous wall is 


V 


s 


dcp 



(53) 


Since the fact that <p is constant along the upper boundary implies that d(p/dX - 0, 


dW 


_i ilk + if] 

i d£ + 3g 

d(p 

dZ 

L 

ro 

ro 

ax dY_ 

dY 



<p=i 


<P= 1 


It therefore follows from equation (50) that 


y = V 

* (Po - Ps> S S 


1 /e^ + A 2 \ 1/2 
A \ e ^ + 1 / 


(54) 


Solution for a Step in Wall Temperature 

In order to demonstrate how equation (42) can be applied, specific results will now 
be obtained for the example of a step function wall temperature variation along the 
boundary S. This solution also yields for zero step height the case of uniform wall 
temperature. In particular, let t g change from t^ to t 2 at an arbitrary location L g 
so that the boundary condition (14a) becomes 
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(55) 


T = 1 + NF(L S - Lg) 

where F is the unit step function defined by 

r F(X) = 1 X > 0 

< 

F(X) = 0 X<0 

The correspondence between L and ii/ and L' and \I/ is obtained from equations (51) 

s s 

and (52). However, because of the special properties of the step function it is not neces- 
sary to relate L g to \f/ in order to evaluate H(/3). In fact, equation (39) becomes in 
this case 


H(jS) = f F ty- !//')e" i/3 ^ dxp 

J-O O 


or by letting X = xj/ - xf/ 


H03) = e 


-x/3^' f F (\)e~ ll3X dX 

J- OO 


The evaluation of this integral is given in reference 7 as 

H(/3) = e" i/3 ^ [Jr6(/3) + 


(56) 


where the notation P. V. indicates that the Cauchy principal value of the integral will 
be taken when integrating l/i/3. 

This expression for H(/3) can now be substituted into equation (42). Upon substituting 
the first term of equation (56) into the integral which appears in equation (42), we find 
that 



(B) Va 2 + e 2 


A cosh(VA 2 + /3 2 ) + Va 2 + (3 2 sinh^V A 2 + /3 2 ) 


A sinh 


inh(f / A 2 ^ + /3 2 )+ >^A 2 + (3 2 cosh(VA 2 + (3 2 


d^A 

2 


(57) 
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Substituting the remaining term of equation (56) into this integral yields the expression 


i/3 


A coshftA 2 + j3 2 )+ r A 2 + /3 2 sinhf f A 2 + /3 2 


A sinhl ? A 2 + /3 2 ) + Va 2 + /3 2 cosh( ?A 2 + 0 2 


As /3 goes to +°° and the integrand approaches — ^ and — e 1 ^^ -1 ^ 

i i 

respectively so that the integral must be interpreted as the Fourier transform of a dis- 
tribution. For this reason this expression is not suitable for numerical evaluation. 
However, in order to obtain an alternative expression for this integral which is suitable 
for this purpose, we need only add and subtract the term 


JL 

2o 


where E (j3) is the sign function which is defined by 


E(j8) = 1 


E03) =-l 


13 >0 
(3 < 0 


It is found from a table of Fourier transforms of distributions that 


\f 


E^eW-^ d/3 




Hence the integral (58) can be written as 
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Upon changing the variable of integration from /3 to -/3 in the first integral, we 
find that the imaginary parts as well as the singularity at the origin cancel out, and we 
obtain 


N 2 + N 

2n (if/ - if/') v 


f 


sin fi(if/ - if/') 


Vx 2 + 1 


X cosh Vx 2 + f? * /x 2 + g 2 sinhVx 2 ♦ 


X sinh 


Vx 2 + /3 2 + Vx 2 + / 3 2 cosh Vx 2 + /3 2 _ 


-l>dj8 


or 


N + N 


n(xf/ - if/) it 


f 


(X - y3) Vx 2 + /3 2 cosh Vx 2 + /3 2 + (X 2 + /3 2 - X/3)sinhVx 2 + /3 2 
X sinh Vx 2 + (3 2 + Vx 2 + /3 2 cosh Vx 2 + /3 2 


X 


sin[/3(y/ - I//')] d/3 {59) 

P 


The integral in this expression is absolutely convergent, and hence, can easily be eval- 
uated numerically. The results contained in relations (57), (59), and (50) show that 
equation (42) yields the following expression for the heat flux entering the boundary S: 


% h r 


1 (e 71 * + A 2 


1/2 


k m (t l _t =o) A V^ 


2X + NX 


e * + 1 


" tF < L s - L s> 


N 


n(if/ - if/') 


+ 

7T 


f 

•4* 


N / (X - (3) Vx 2 + /3 2 cosh Vx 2 + /3 2 + (X 2 + /3 2 - X/3)sinhVx 2 + /3 2 


,2 n 2 1 . ^ ^2 n 2 . / x 2 n 2 

X sinh Vx 2 + 0 2 + }/ )? + (3 2 coshVx 2 + /3 2 


sin[/3(^/ - i//)] d/3 l 

P 


(60) 


25 


I 


Solution for Uniform Heat Flux 


In order to illustrate how equation (45) is applied, consider the case where a uniform 
heat flux is imposed along the surface S. The function G in equation (6b) is zero in this 
instance, and equation (50) gives an expression for dW/dZ . Hence in this case the 

function E in equation (44) is given by 



As \p — to the integrand goes to e and as ip — -«> the integrand goes to 
e^/A, so that this integral exists only as the Fourier transform of a distribution and 
not as an ordinary integral. In order to separate out the singular parts of this integral, 
notice that equation (61) can be written as 



where the first two integrals are absolutely convergent and the last two integrals can be 
evaluated from a table of Fourier transforms to give 

1 2*5 03) + ^1 - -i^|*6G3) + 


Thus the following expression for E (j3) is obtained 
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This is now inserted into equation (45). We shall consider the integrals resulting 
from the various terms of equation (62) individually. The integral resulting from the first 
term of equation (62) is 


A 

2 


H-i)f 

CO 


6 03 ) 


A sinh^A 2 + /3 2 ) + V^A 2 + /3 2 cosh^A 2 + j3 2 ) 


(2 A 2 + ;3 2 )sinh( }/ A 2 + /3 2 ) + 2A^A 2 + /3 2 cosh(^A 2_ T/3 2 / 

X sinh A + A cosh X A + 1 


(fi) 


2A 2 sinh A + 2A 2 cosh A 


4A 


e i/3 ^ d/3 


(63) 


The integral resulting from the second term of equation (62) is 


(A - 1) 

2n 


P. V. 


f 


A sinh^A 2 + 0 2 ) + Va 2 + /3 2 cosh( V A 2 + /3 2 


>i/3^ 


(2 A 2 + /3 2 )sinh^ V A 2 + /3 2 ) + 2A V >- 2 + /3 2 cos!i(Va 2 + /3 2 ) ^ 


A + j3 coshy r A + /3‘ 


d/3 


= (A - 1) 

IT 


f 


sinh^V A 2 + j3 2 ) + Va 2 + /3 2 cos1i(Va 2 + /3 2 ) 


(2 A 2 + /3 2 )sinh( Va 2 + /3 2 ) + 2aVa 2 + ° 2 — •-f ^' 2 ■ ° 2 




0 coshl f A + f3‘ 


sin^d/3 


I 3 


(64) 


This latter integral is absolutely convergent. Now consider the integrals resulting from 
inserting the last two terms of equation (62) into equation (45). These integrals are also 
found to be absolutely convergent. The exponential terms in the integrals can therefore 
be expressed in trigonometric form. After collecting real and imaginary parts, it is 
found that the imaginary part of the integrand is an odd function of /3 so that its integral 
from /3 = -<=o to °° is zero. The integral of the remaining real part can be expressed as 
an integral from /3 = 0 to °° and is then combined with equations (63) and (64) to yield the 
following expression for the surface temperature: 


(t s too ^ k m A + 1 


4A 


«s h r 


V 


E(/3, 4s) 


V'a 2 + /3 2 cosh ^A 2 + /3 2 + A sinh 


inhV^A 2 + B 2 


2A 




/3 2 cosh }/ A 2 + /3 2 + (2 A 2 + /3 2 )sinhVA 2 + /3 2 


d/3 (65a) 
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where 



DISCUSSION 

A general analytical method has been presented for determining the heat transfer 
characteristics of two-dimensional porous cooled media. This method is based on trans- 
forming the energy equation so that its independent variables are the coordinates in the 
potential plane for the fluid flow. The energy equation is separable in these coordinates, 
and hence a general solution can be obtained. A feature that aids in the solution is that 
the boundaries of the porous material are at constant pressure and as a consequence they 
map into parallel constant velocity potential lines in the complex potential plane. This 
provides a convenient region in which to solve the energy equation. The solution in the 
potential plane is then related to the physical plane by conformal mapping between the two 
regions. 

The following discussion concerns the application of the solution to the example of a 
porous wall with a step change in cross section. No attempt is made to conduct a para- 
metric study to demonstrate the effect of the governing dimensionless parameters since 
this would only be of practical value if a particular engineering application were in mind. 

Consider now the solution given by equations (51), (52), and (60). They give the 
dimensionless heat flux into the boundary S of the porous material parametrically as a 
function of position when the temperature distribution is a step function along that sur- 
face. The result contains the parameters X, N, i//, and A. The parameter A is the 
ratio of the thicknesses of the two regions of the porous wall as shown in figure 4(a). The 
only results which will be presented are for a value of A equal to 2. Hence, the geom- 
etry will always be as shown in figure 6(a). The heat flux is first computed as a function 
of if/ from equation (60). It is then expressed as a function of the coordinate L along 
the surface by use of the correspondence between L g and if/ given by equations (51) 
and (52). 

The dimensionless velocity leaving the upper surface is given by equation (54) and is 
shown in figure 6(b) for the porous plate with a step ratio A of 2. At the limits of large 
negative L g and large positive L g the velocity V g goes, respectively, to 1 and 1/2, as 
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(a) Geometry of porous medium. 



(b) Dimensionless exit velocity from upper surface of medium. 



Figure 6. - Heat transfer to porous medium with upper 
surface at constant temperature, A * 2, N = 0. 


would be expected by examining the one- dimensional solution given in reference 1. The 
extent of the two-dimensional region required for the velocity to change between the two 
limiting values is shown by the figure. 

If the parameter N = (t 2 - tj)/^ - t^) is set equal to zero, then t 2 = tj and the 
surface S is at uniform temperature. With N = 0 the solution in equation (60) becomes 

% h r _ 1 (e^K A 2 \ 1/2 2x 
k m (t l _ to °) A \ e ^ + 1/ 


so it follows from equation (54) that 


% h r 


2 Vi-y 
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Thus for a given porous matrix when the surface temperature is uniform, the dimension- 
less heat flux into the boundary is directly proportional to the fluid exit velocity. The 
dimensionless heat flux is shown in figure 6(c). This curve is the same as the curve in 
figure 6(b). 

The results of the calculations for N = 1 and A = 1 are shown in figure 7. In this 
case the surface temperature has the distribution shown in figure 7(b). The heat flux is 
computed for three different positions of the step corresponding to L* = -1, 0, +1. The 

b 

values of i//' which are the values of corresponding to these values of L are found 

s 

from equations (51) and (52), and then used in the solution equation (60). For large L 

s 

the dimensionless heat flux is found to be the same as for large negative L_, as shown 

s 



(a) Geometry of porous medium. 



0 ' — 

(b) Temperature distributions along upper surface. 



(c) Dimensionless heat flux at upper surface correspond- 
ding to temperature distributions in part (b). 

Figure 7. - Heat transfer to surface for various locations 
of step in surface temperature, A = 2, A = 1, N = 1. 
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by figure 7(c). This is the result of the compensating effects of having twice the temper- 
ature drop and half the exit velocity in the thick region of the wall as compared with the 
thin region. In the general case these two limiting values of q g are not equal. In the 
vicinity of the temperature jump q changes discontinuously from a large negative value 
to a large positive value . This is a local heat conduction effect wherein the heat which is 
conducted into the porous wall at the high temperature side of the step flows into the low 
temperature region, where it is then extracted from the wall to maintain the surface tem- 


perature discontinuity. 

Figure 8 shows the effect of varying the parameter A = — — 

2k m 


while the 



(a) Geometry of porous medium. 


2 1 — 

i f J 

urt i — t 

0 — 

(b) Temperature distribution along upper surface. 



L $ = l s /h r 


(c) Dimensionless heat flux at upper surface. 

Figure 8. - Effect of X on heat transfer to surface of 
step porous medium with step in surface temperature, 
A - 2, N - 1. 
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temperature distribution on the wall surface remains fixed as given in figure 8(b). The 
larger values of X are associated for example with larger pressure differences p - p 
and hence with larger flows. Thus to maintain a fixed surface temperature, the heat flux 
into the plate can be increased when the flow is increased to correspond to a larger value 
of X. 

Now consider the solution given by equations (51), (52), and (65) for the case of uni- 
form heat input at the surface S. There are only two parameters involved, A and X. 
Numerical results have been obtained for A = 2, and they are shown in figure 9. The 

exit velocity V_ is the same as in figure 6 and is repeated here for convenience in in- 
s 



(a) Geometry of porous medium. 



(c) Temperature of upper surface. 

Figure 9. - Effect of A on surface temperature for im- 
posed uniform heat flux, A = 2. 
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terpreting the results. Since X increases when the flow through the wall increases, the 
temperature difference t - t^ between the surface and the reservoir decreases as X 
is increased for a fixed q_. Note that q_ has been combined into the dimensionless 
ordinate of figure 9(c) so that each curve can be used for any value of q„. As would be 

o 

expected, the regions of high surface temperature are associated with the regions of low 
exit velocity. 

The preceding results serve to illustrate the type of heat transfer behavior to be ex- 
pected in a two-dimensional porous configuration. 


CONCLUSIONS 

An analytical approach has been developed for obtaining the heat transfer behavior 
in a two-dimensional porous material. The method depends on the fact that the surfaces 
of the porous material at the inlet and outlet of the coolant are each at constant pressure. 
For the type of flow involved, the dimensionless pressure can be regarded as a velocity 
potential, and as a consequence the porous region maps into a simple strip in the com- 
plex potential plane. The energy equation was transformed into the potential plane, and 
a general solution was obtained. This solution can be applied to any wall geometry by 
finding the appropriate conformal map of the strip in the potential plane into the physical 
plane . 

The analytical technique was applied for various imposed thermal conditions in the 
specific case of a step porous wall which is made up of two regions, each having a differ- 
ent uniform thickness. 
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